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Abstract: Presented are complete next-to- leading order electroweak (NLO EW) correc- 
tions to top-squark pair production at the Large Hadron Collider (LHC) within the Minimal 
Supersymmetric Standard Model (MSSM). At this order, also effects from the interference 
of EW and QCD contributions have to be taken into account. Moreover, photon-induced 
top-squark production is considered as an additional partonic channel, which arises from 
the non-zero photon density in the proton. 
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1. Introduction 

Within supersymmetric theories top-squarks are the supersymmetric partners of the left- 
and right-handed top quarks. The two superpartners ti and tn, which belong to chiral 
supermultiplets Q and T, in general mix to produce two mass eigenstates t± and £2- hi 
many supersymmetric models the lighter mass eigenstate appears as the lightest colored 
particle for reasons related to the large top Yukawa coupling. The large mixing in 
the stop sector leads to a substantial splitting between the two mass eigenstates, and the 
evolution from the GUT scale to the electroweak scale yields low values for the stop masses 
when a universal scalar mass is assumed at the high scale The search for top-squarks 
is therefore of particular interest for the coming LHC experiments, where they would be 
primarily produced in pairs via the strong interaction, with relatively large cross sections. 

Current experimental limits on top-squark pair production include searches performed 
at LEP HI reviewed e.g in Q, and at the Tevatron, done by the CDF and D0 collaborations 
in approximately 90 pb _1 of Run I data |J. Extended searches have been done using Run II 
data samples by both CDF and D0 ||. Limits on the top-squark mass, depending on the 
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mass of the lightest neutralino, are provided with the assumption that BR(i\ — ► cxi) = 
100 % in §. 

Experimental searches for the top-squarks have also been done in ep collisions at 
HERA (8|, where only single stop production could be kinematically accessed and hence 
constraints have been derived essentially on the R-parity violating class of supersymmetric 
models. 

Concerning the theoretical predictions, QCD-based Born-level cross sections for the 
production of squarks and gluinos in hadron collisions have been calculated in ||. They 
have been improved by including NLO corrections in supersymmetric QCD (SUSY-QCD), 
worked out in with the restriction to final state squarks of the first two generations, 
and for the stop sector in [jllj] . The production of top-squark pairs in hadronic collisions is 
diagonal at lowest order at 0(a 2 ). Electroweak (EW) contributions of (D(a 2 ) are suppressed 
by two orders of magnitude. Also at 0(o^) the production mechanism is still diagonal. 
Non-diagonal production occurs at 0(ctg), and the cross section is accordingly suppressed. 
Production of non-diagonal top-squark pairs can also proceed at 0(a 2 ) mediated by Z- 



exchange through qq annihilation [|12| as well as in e + e annihilation |13 |. 

The LO cross section for diagonal top-squark pair production depends only on the mass 
of the produced squarks. As a consequence, bounds on the production cross section can 
easily be translated into lower bounds on the lightest stop mass. At NLO, the cross section 
becomes considerably changed and dependent on other supersymmetric parameters, like 
mixing angles, gluino mass, masses of other squarks, etc., which enter through the higher 
order terms. Once top-squarks are discovered, measurement of their masses and cross 
sections will provide important observables for testing and constraining the supersymmetric 
model. 

In the following, we study the NLO contributions to diagonal top-squark pair pro- 
duction that arise from the electroweak interaction within the Minimal Supersymmetric 
Standard Model (MSSM). We assume the MSSM with real parameters, R-parity conser- 
vation, and minimal flavor violation. The outline of our paper is as follows. In Section 2, 
we present analytical expressions for the partonic and hadronic LO cross sections. We also 
introduce some basic notations used throughout the paper. Section 3 is dedicated to the 
classification of the NLO EW contributions into virtual and real corrections with the treat- 
ment of soft and collinear singularities, and photon-induced contributions. In Section 4, 
we give a list of input parameters and conventions, followed by our numerical results for 
the hadronic cross sections and distributions for pp collisions at a center-of-mass energy 
■*/S = 14 TeV at the LHC. We also investigate the application of kinematical cuts, and we 
analyze the impact of varying the MSSM parameters. 



2. Top-squark eigenstates and LO cross sections 

In the MSSM Lagrangian, mixing of the left- and right-handed top-squark eigenstates t^m 
into mass eigenstates t\ji is induced by the trilinear Higgs-stop-stop coupling term A t and 
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the Higgs- mixing parameter [i. The top-squark mass matrix squared is given by [14] 

an- - ( '"' „ ! . (2.1) 



.2 ( m t + a ll m t B LR \ 
\m t B LR mf + Cijfl/ ' 



with m( denoting the top-quark mass and 

1 



A LL = ( 2 ~ 3 sin ® w ) m ^ cos 2 ^ + m 
-Blh = A t - (i cot /3 , 

Crr = o sm2 m z cos 2/3 + f^ 2 - 



(2.2) 



3 ' ' ' 

Here, tan (3 is the ratio of the vacuum expectation values of the two Higgs doublets and 
m <5 3 > m u 3 are ^ ne soft-breaking mass terms for left- and right-handed top-squarks, respec- 
tively. 

The top-squark mass eigenvalues are obtained by diagonalizing the mass matrix, 

,2 
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and the mixing angle Or is determined by 



tan 2% 
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(2.3) 
(2.4) 

(2.5) 



All — Crr 

At hadron colliders, diagonal pairs of top-squarks can be produced at leading order in 
QCD in two classes of partonic subprocesses, 



(2.6) 



99 -> h i\ and t 2 t* 2 , 
qq — ► ti t\ and t 2 i* 2 , 

where qq denotes representatively the contributing quark flavors. The corresponding Feyn- 



man diagrams for the example of t\ t* production are shown in the appendix, Fig. |A.l 



As already mentioned, mixed pairs cannot be produced at lowest order since the g tt* and 
ggtt* vertices are diagonal in the chiral as well as in the mass basis. 
The differential partonic cross sections for the subprocesses, 

1 



da 99m (s) 



J2\M 9 9 ' q \s,t,u)fdt 



(2.7) 
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can be expressed in terms of the squared and spin-averaged lowest-order matrix elements, 
as explicitly given by [pXl], 
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with t r = t — m~ , u r = u — mf , where s, i, -u are the usual Mandelstam variables, i = 1, 2 
denotes the two mass eigenstates. The 577(3) color factors are given by A^ = 3, Co = 
AT(At2 _ i) = 24, C x = (N 2 - 1)/N = 8/3 and C F = (N 2 - 1)/(2N) = 4/3. 

The differential cross section at the hadronic level for the process AB — * tit*, i =1, 2, 
is related to the partonic cross sections through 

rl jrAB 

da AB (S) =J2 dr^- da?(s) , (2.10) 

a,b Jr ° aT 

with r = s/S, S (s) being the hadronic (partonic) center-of-mass energy squared and 
ro = 4m? / S is the production threshold. The sum over a, b runs over all possible initial 
partons. The parton luminosities are given by 



dC#*_ 1 f x dx 



dr 1 + 



J ^ fa/A^^F) fb/B^VF) + fb/A^* VF*) /o/B^j^f) 



(2.11) 



where the parton distribution functions (PDFs) f a /A( x i parameterize the probability 
of finding a parton a inside a hadron ^4 with fraction x of the hadron momentum at a 
factorization scale \ip. 

3. Classification of EW NLO corrections 

In the following we describe the calculation of EW contributions to top-squark pair pro- 
duction at NLO. For the treatment of the Feynman diagrams and corresponding ampli- 
tudes we make use of FeynArts 3.2 [1^] and FormCalc 5.2 with LoopTools 2.2 | |16| |, 
based on Passarino-Veltman reduction techniques for the tensor loop integrals which 



were further developed for 4-point integrals in IS]. Higgs properties are computed with 



FeynHiggs 2.5.1 [19|. 



The super symmetric final state does not allow to separate the SM-like corrections from 
the superpartner contributions which are necessary for the cancellation of ultra-violet (UV) 
singularities. As the photino is not a mass eigenstate of the theory, it is also not possible 
to split the EW corrections into a QED and a weak part, which is often the case in SM 
processes. In order to obtain a UV finite result, we have to deal with the complete set of 
EW virtual corrections including photonic contributions. These are infrared (IR) singular 
and thus also the real photonic corrections have to be taken into account. In addition, 
a photon-induced subclass of corrections appears at NLO as an independent production 
channel. 

3.1 Virtual corrections 

The virtual corrections arise from self-energy, vertex, box, and counter-term diagrams. 
These are shown in the appendix, in Fig. A. 3 for the qq annihilation and and in Fig. A. 4 for 



the gluon fusion channel, respectively. Getting an UV finite result requires renormalization 
of the involved quarks and top-squarks. The renormalized quark and squark self-energies 
are obtained from the unrenormalized initial quark self-energies 

= ^l(p 2 ) +^ + ^ R (p 2 ) + m q m P 2 ), (3.1) 
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according to 



±1(P 2 ) = ^ L (P 2 )+6ZI, 

^ 2 ) = S|(p 2 )-i(^l + ^) + ^i 



(3.2) 



and from the top-squark self-energies £ t ~.(/c 2 ) (for i =1, 2), according to 



t k (k 2 ) = S fi (fc 2 ) + k 2 5Z~ u - m\5Z~ ti - 6m\ , 



(3.3) 



with the renormalized quantities denoted by the symbol £. 

The full set of virtual contributions is UV finite after including the proper counter- 
terms for self-energies, quark vertices, and squark triple and quartic vertices, as listed in 
the following set of Feynman rules: 



V 

ti 



z<5£ t ~ = i (k 2 8Z~ ti - m\SZ- u - <5m 2 - ) , (3.4) 
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^ = -ig s T c {k + k') ti 5Z k , 



(3.5) 
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1 igl Q<U + d ahc T c ^j 9ilv 5Z~ ti , (3.6) 



iSAl = -ig s T c ^ (uSZl + u+6Z R ) , (3.7) 



where k, k! denote the momenta of top-squarks (in the direction of arrows), a, b, and 
c are the gluonic color indices, T c and d a i, c are the color factors (we skip the fermionic 
and sfermionic color indices), and u>± = (1 ± 7s)/2 are the projection operators. The 
renormalization constants are fixed within the on-shell renormalization scheme as follows, 



5m-, 
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(3.8) 
(3.9) 

(3.10) 
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There is no renormalization of the gluon field at 0(a). Also, the strong coupling 
constant does not need renormalization since UV singularities cancel in the sum of 3- 
and 4-point functions and their corresponding counter-terms from quark and squark field 



renormalization (see Figs. A. 3 and A. 4 in the appendix). 

Loop diagrams involving virtual photons generate IR singularities. According to Bloch- 
Nordsieck @, IR singular terms cancel against their counterparts in the real photon 



corrections. To regularize the IR singularities we introduce a fictitious photon mass A. In 
case of external light quarks, also collinear singularities occur if a photon is radiated off a 
massless quark in the collinear limit. We therefore keep non-zero initial-state quark masses 
m q in the loop integrals. This gives rise to single and double logarithmic contributions of 
quark masses. The double logarithms cancel in the sum of virtual and real corrections, 
single logarithms, however, survive and have to be treated by means of the factorization. 

In the gg fusion channel, IR singularities originate only from final-state photon ra- 
diation, and mass singularities do not occur. In the qq annihilation subprocess, the IR 
singular structure is extended by the contributions related to the gluons which appear in 
the 4-point UV finite loop integrals. There are two types of IR singular box contributions 
(Fig. A.3| c). The first group is formed by the gluon-photon box diagrams with two sources 



of IR singularities, one related to photons, the other to gluons. The second group consists 
of the gluon- Z box diagrams with IR singularities originating from the gluons only. There 
is also an IR finite group of 0(aa s ) box diagrams which consists of gluino-neutralino loops 



(Fig. A.3| d). Owing to the photon-like appearance of the gluon in the box contributions, 



the gluonic IR singularities can be handled in analogy to the photon IR singularities. 
3.2 Real corrections 

To compensate IR singularities in the virtual EW corrections, contributions with real 



photon (Fig. A.5 a and c) and real gluon radiation are required. In case of gg fusion, 
only photon bremsstrahlung is needed, whereas in the qq annihilation channel, also gluon 
bremsstrahlung at the appropriate order 0(aa 2 ) has to be taken into account (Fig. A. 6) 



to cancel the IR singularities related to the gluon. The necessary contributions originate 
from the interference of QCD and EW Born level diagrams, which vanishes at LO. Not 
all of the interference terms contribute. Due to the color structure, only the interference 
between initial and final state gluon radiation is non-zero. 

Including the EW-QCD interference in the real corrections does not yet lead to an IR 
finite result. Also the IR singular QCD-mediated box corrections interfering with the 0(a) 
photon and Z-boson tree-level diagrams are needed. Besides the gluonic corrections there 
are also the IR finite QCD-mediated box corrections, which contain gluinos in the loop. 
Interfered with the 0{a) tree-level diagrams, these also give contributions of the respective 
order of Oiaaz ). The set of all 0{a 2 s ) dia grams is shown in Fig. [A.7|. 



So far we have mentioned only the IR singular bremsstrahlung contributions. However, 
there are also IR finite real corrections to both gluon fusion and qq annihilation processes. 
In addition to the photon radiation off the off-shell top-squark there are photon radiation 
contributions originating from the quartic gluon-photon-squark-squark coupling. These 



contributions do not have to be regularized since they are not singular (Fig. A.5| b and d) 
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The treatment of IR singular bremsstrahlung is done using the phase space slicing 
method. Imposing cut-offs AE on the photon/gluon energy and A9 on the angle between 
the photon/gluon and radiating fermion, the photonic/gluonic phase space is split into soft 
and collinear parts which contain singularities and a non-collinear, hard part which is free 
of singularities and is integrated numerically. The sum of virtual and real contributions, 
each of them dependent on the cut-off parameters AE and A9, has to provide a fully 
independent result. To ensure this we perform numerical checks. 

In the singular regions, the squared matrix elements for the radiative processes fac- 
torize into the lowest-order squared matrix elements and universal factors containing the 
singularities. 



3.2.1 Soft singularities 

The soft-photon part of the radiative cross section in the qq annihilation channel 

<V '4 Coft + e? 5{T ft + 2e q e t S%%) d&f(§) , (3.11) 



7T 



and in the gg fusion channel 
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(3.12) 



can be expressed using universal factors, S so ^ ' , which refer to the initial state radiation, 
final state radiation or interference of initial and final state radiation, respectively, da^ 99 
denote the corresponding partonic lowest order cross sections. The singular universal fac- 
tors, similar to those in [Ell, read as follows, 
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(3.13) 



Here, e q and ej are the electric charges of the initial quark and of the top-squark, re- 
spectively, and we introduced S s = 2AE/\/H, where AE is the slicing parameter for the 
maximum energy a soft photon may have. For application purposes, it is useful to express 
Eq. (3.13) in terms of Mandelstam invariants, t and u, using the relations 



t, u 



(1 =F Pcost 
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(3.14) 
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The soft-gluon part for the qq channel can be written in a way similar to ( |3.11| ), but 
with a different arrangement of the color matrices, 



d&llft,g(s) — — fisoft ^TfjTjiT^Tml 



7T 

X 



2ReJ2(MZ*M™+M™*Ml% 



dt 



(3.15) 



16tts z 



with A4 denoting the "Born" matrix elements for g, 7 and Z exchange where the color 
matrices are factorized off. Explicitly, it can be written as follows, 



8e q e t (JJJii) 2 -2e t sm 2 6 w )(e-4e q sm 2 6 w ) 



16ir 2 aa s 
4-9 



(i-m k )(u-m~ k ) 



sin 2 9w cos 2 6w s(s — m 2 z ) 
di 



(3.16) 



2 \ 2 ~ 
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involving the top-squark mixing matrix of Eq. ( |2,4| ), and e = ±1 for up- and down- type 
initial quarks, respectively. 



3.2.2 Collinear singularities 

Collinear singularities arise only from initial-state photon radiation in qq annihilation. The 
collinear part of the 2^3 cross section is proportional to the Born cross section of the 
hard process with reduced momentum of one of the partons. Assuming that parton a with 
momentum p a radiates off a photon with p 7 = (1 — z)p a , the parton momentum available 
for the hard process is reduced to zp a . Accordingly, the partonic energy of the total process 
inclusive photon radiation is s = (p a + Pb) 2 = fS , and for the hard process the reduced 
partonic energy is s = (zp a + p^) 2 = tS . The 'total' and 'hard' variables are thus related 
by s = zs and t = zf . 

Having defined these variables, the partonic cross section in the collinear cones can be 
written in the following way [22, 23] 



dcr co u{s) 
with Kcoll(z) 
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(3.17) 



where P qq {z) = (1+z )/(l — z) is an Altarelli-Parisi splitting function [24] and 5g is the cut- 



off parameter to define the collinear region by cos 6 > 1 — 5g. The Born cross section refers 
to the hard scale s, whereas in the collinear factor the total energy s is the scale needed. 
In order to avoid an overlap with the soft region, the upper limit of the z-integration in 
Eq. (3.17) is reduced from z = 1 to z = 1 — 8 S . 

As already mentioned, after adding virtual and real corrections, the mass singularity in 
Eq. ( 3.17] ) does not cancel and has to be absorbed into the (anti-)quark density functions. 
This can be formally achieved by a redefinition of the parton density functions (PDFs) at 



- 8 - 



NLO QED as follows [||, ||, [26 
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The QED factorization scheme dependent A sc -parameter is X sc = in the MS'-scheme and 
A sc = 1 in the DIS scheme. 

At the hadronic level, we define the collinear part of the real corrections for the case 



where parton a radiates off a collinear photon, in the following way by use of Eq. (3.18|), 
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(3.19) 



where the lower limit of the z-integration is constrained to x, since the parton momentum 
fraction x/z has to be smaller than unity. The integral is free of any mass singularity, 
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(3.20) 



The K^Jjf-teim in Eq. ( |3,lSj ) cancels the mass singularities owing to soft photons that 
remain in the sum of the virtual corrections and the soft correction factor 5™j t in Eq. ( 3.13j ) . 



3.3 Photon-induced top-squark pair production 

We also consider the photon-induced mechanisms of the top-squark pair production. At 
the hadronic level, these processes vanish at leading order owing to the non-existence of a 
photon distribution inside the proton. At NLO in QED, however, a non-zero photon density 
arises in the proton as a direct consequence of including higher order QED effects into the 
evolution of PDFs, leading thus to non-zero photon-induced hadronic contributions. 

Feynman diagrams corresponding to the photon-gluon partonic process are illustrated 



in Fig. A. 2 . Although these are contributions of different order, they are tree- level contribu- 
tions to the same hadronic final state and thus deserve a closer inspection. The differential 
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cross section for this subprocess is 



da^(s) 
EN 
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(3.21) 



expressed in terms of the reduced Mandelstam variables t r = t — m~ , u r = it — mh . The 
quark-photon partonic processes represent contributions of higher order and we do not 
include them in our discussion here. 

The photon density is part of the PDFs at NLO QED, which have become available 
only recently f^Tj; here we present the first study of these effects on the top-squark pair 
production. 



4. Numerical results 

For the numerical discussion we focus on the production of light top-squark pairs t*ti in 
proton-proton collisions for LHC energies. We present the results in terms of the following 
hadronic observables: the integrated cross section, a, the differential cross section with 
respect to the (photon inclusive) invariant mass of the top-squark pair, (da/dMi nv ), the 
differential cross sections with respect to the transverse momentum, (dcr/dpr), to the 
rapidity, (da/dy), and to the pseudo-rapidity, (da/drj), of one of the final state top-squarks. 
For getting experimentally more realistic results for the cross sections we also apply typical 
sets of kinematical cuts. A study of the dependence on the various SUSY parameters is 
given towards the end of this section. 

The NLO differential cross section at the hadron level is combined from the contribut- 
ing partonic cross sections by convolution and summation as follows, 

r 1 ( dL vv - rir pp rir pp 1 

da pp (S) = J dr I £ da^(s) + ^ da^(s) + ^ da 9 ^s)\ , (4.1) 

where da qi ^ and da" represent full one-loop results, including complete virtual and real 
corrections, and da 9 ^ is given in Eq. ( 3.21[ ). The respective parton luminosities refer 



to Eq. (pUD . 

One has to take care of the fact that each top-squark observed in the laboratory 
system under a certain angle 9 can originate from two different constellations at parton 
level: parton a(b) out of hadron A(B) and vice- versa, corresponding to 9 — > (ir — 9). 
Both parton level configurations have to be added correctly for hadronic distributions (for 
explicit formulas see e. g. p8| ). Note that the two boost factors j3 relating the two partonic 
center-of-mass (c. m. ) systems with the laboratory system differ by a relative sign, as do 
the rapidity and the pseudo-rapidity of each particle. 

Assuming that the forward-scattered parton a carries the momentum fraction x of 
hadron A and the backward-scattered parton b the momentum fraction t/x of hadron B, 
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the boost factor (3 is given by 



P= X —^. (4.2) 

X + T/X 



The rapidity of one of the final state top-squarks in the laboratory system, y(= y$*), is 
related to the rapidity in the partonic cm. frame, y cm = artanh(p^ m / E cm ) , via a Lorentz 
transformation, 

y = y cm - artanh(-/3) = y cm + - In — . (4.3) 

2 r 

The pseudo-rapidity rj is related to r] cm = — ln(tan 9 cm /2) in the cm. frame via 



„ = arsinh | ^<A + cosl S, fm (^-^) + isinh,™ (^ + ^| | , (4.4) 
which can be derived using the representation 



p = y\J m i +PyCOsh t/, 0, pr> PTsinh^J (4.5) 

for the top-squark momentum p = p^*. Since the final state particles are massive, rapidity 
and pseudo-rapidity do not coincide; in the limit m — * one obtains rj = y. 

4.1 Input parameters and conventions 

Our Standard Model input parameters are chosen in correspondance with fEfl , 

M z = 91.1876 GeV, M w = 80.403 GeV, 
a' 1 = 137.036, a{M z y x = 127.934, G F = 1.1664 x 10" 5 GeV" 2 , (4.6) 
m t = 172.7 GeV, m b = 4.7 GeV, m b (m b ) = 4.2 GeV . 

All lepton and all other quark masses are set to zero unless where they are used for regu- 



larization. As a reference we consider the SPA SUSY parameter point SPS la' [29|, unless 
stated otherwise. The current value of the top-quark mass, mt = 170.9 ± 1.9 GeV p0| , 
increases the top-squark mass by 0.2%, which reduces the total cross section by w 1%. 
The changes for the relative corrections are completely negligible. 

For the parton distributions, we use the set MRST 2004 QED |27]], as already mentioned 
previously. Factorization and renormalization scales are chosen equal, p,p = [ir = 2m^. 

4.2 Hadronic cross sections and distributions 

In Table |l] we show results for the cross section for top-squark pair production at the 
LHC within four different scenarios, chosen out of the SPS benchmark scenarios of the 



minimal SUGRA type [29, 31 1. The integrated hadronic cross sections at leading order, 
a LO , the absolute size of the EW corrections corresponding to the difference between the 
LO and NLO cross sections, Aa NLO , and the relative corrections, 5, given as the ratio of 
NLO corrections to the respective LO contributions, are presented for the gg fusion, the 
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scenario 


channel 


a 


LO rrui 
|lbj 


A NLO T-fUl 

Act ^ [ibj 


r Aa NLO 




SPS la 
= 376.2 GeV) 


qq 
99 
91 
total 


222 
1444 

1666 


(+0.985) 


-9.71 
-15.4 
29.0 
3.90 


-4.4% 
-1.1% 

023% 




SPS la' 

= 322.1 GeV) 


qq 
99 
91 
total 


439 
3292 

3731 


(+1.88) 


-11.6 
-14.6 
58.5 
32.3 


-2.6% 
-0.44% 

0.87% 


( m tl 


SPS 2 
= 1005.7 GeV) 


qq 
99 
91 
total 


1.17 
2.97 

A 1 A 

4.14 


(+0.00539) 


-8.99 x 10~ 2 
-3.07 x 10~ 2 
15.5 x 10~ 2 

o.44 X 1U 


-7.7% 
-1.0% 

0.8o7o 


( m h 


SPS 5 
= 203.8 GeV) 


qq 
99 
91 
total 


2900 
31960 

34860 


(+10.2) 


-13.3 
499 
405 
891 


-0.46% 
1.6% 

2.6% 



Table 1: Numerical results for the integrated cross sections for light top-squark pair production 



at the LHC within different SPS scenarios |29, 31 



qq annihilation, and the #7 fusion channel separately. The 37 channel contributes only 
at NLO. For the qq channel, also the numbers for the 0(a 2 ) pure electroweak Born level 
contributions are given in brackets. These are typically smaller by one order of magnitude 
compared to the EW NLO corrections. 

In scenarios where the top-squark ti is of intermediate or high mass (as SPS la, SPS la', 
and SPS 2) the NLO contributions are below 1%. The corrections to the qq and the gg 
channels are negative, whereas the 37 contribution is always positive and of the same size 
as the other corrections or even larger. The situation is different in scenarios where the 
top-squark is very light, i.e. lighter than half of m^o, the mass of the heavier neutral Higgs 
boson H°, where a large fraction of the squarks appears through production and decay of 
H° particles. This is the case in the SPS 5 scenario [m { - = 204 GeV, m H o = 694 GeV and 
T(H°) = 9.7GeV derived from FeynHiggs [19|]. The electroweak contributions in the gg 
channel are positive and slightly larger than the 57 fusion contribution. 

The interplay of the three production channels is illustrated in Fig. [l| where the absolute 
EW contributions Aa per channel are shown as distributions with respect to px, Mi nv , y, or 
■q. Owing to the alternating signs, compensations occur where in particular the 57 channel 
plays an important role. 

For realistic experimental analyses, cuts on the kinematically allowed phase space of 
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Figure 1: Comparison of EW NLO contributions from the various parton channels, for the dis- 
tributions of transverse momentum PT{t*), invariant mass of the stop pair, rapidity y(t*), and 
pseudo-rapidity r)(t*) (from upper left to lower right), y and rj are given in the laboratory frame. 
For gg fusion and qq annihilation, A denotes the difference between NLO and LO distributions 
(Act ee Aa NLO ), for gl one has Act ee erg 7 . 

the top-squarks have to be applied. They can be realized by a lower cut on the transverse 
momenta of the final-state particles to focus on high-py jets. Moreover, detectability of 
the final state particles requires a minimal angle between the particles and the beam axis. 
Therefore, we set a cut on the pseudo-rapidity of the top-squarks restricting the scattering 
angle 9 to a central region. Two exemplary sets of cuts are applied in the following figures 
(Figs. 2-5), 

cuts 1: pt > 150 GeV and \rf\ < 2.5 (i. e. 9.4° < 9 < 170.6°), 
cuts 2: p T > 250 GeV and \q\ < 2.5 . 

The differential cross sections and the influence of cuts are the content of Figs. [2] and |3|. 
Displayed are the hadronic cross sections at NLO, differential with respect to pr, Mi nv and 
to y, 77, respectively. Both the full (unconstrained) distributions and the distributions 
with cuts applied are shown. The reduction of the integrated cross section owing to the 
application of cuts is summarized in Table ||. 



-13- 




P T [ GeV l M inv [GeV] 



Figure 2: Comparison of EW NLO differential hadronic cross sections (solid lines) and the dis- 
tributions where kinematical cuts on the final top-squarks are applied for all three production 
channels, gg fusion (upper red plots), qq channels (middle blue plots), and 37 fusion (lower green 
plots). Cuts 1 (dashed lines): p T > 150 GeV, \r)\ < 2.5, cuts 2 (dotted lines): p T > 250 GeV, 
\r)\ < 2.5. Distributions with respect to the transverse momentum px(ti) (left) and the invariant 
mass of the stop pair (right) are shown for t\t\ pair production at the LHC within the SPS la' 
scenario. 

The application of cuts reduces the gg and 37 channels strongly, cutting off the peak 
of the pr-distributions. The reduction is less pronounced in the qq channels where the 
pr-distribution is harder. The px-cuts also shift the threshold of the invariant mass distri- 
butions towards higher values affecting again mainly the gg and 57 channels in height and 
shape. The situation for the rapidity distribution is similar. In the qq channel, the harder 
pr-distribution goes along with a narrower ^-distribution, as shown in the right panels of 
Fig. ||. Most of the top-squarks produced via qq annihilation can be found in the central 
region. In contrast, top-squarks from gg or gj fusion are often produced in the strong 
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Figure 3: Same as Fig. |[ but with respect to the rapidity y(t*) (left) and the pseudo-rapidity 
r,{tl) (right). 

forward (or backward) direction, and the application of cuts on the pseudo-rapidity thus 
reduces the number of gg or g-f based events significantly. 

In order to illustrate the numerical impact of the NLO contributions on the LO cross 
section, we show in Fig. || K factors K = a NLO /a LO for the gg and the qq channel, 
respectively, as distributions with respect to pt and Mi nv . The application of cuts influences 
the K factors only at low values of pt and Mi nv . The EW corrections in the ^-distribution 
reach typically —10% in the gg channel, and —20% in the qq channel, for large values of 
pr- In the invariant mass distributions, they are somewhat smaller, but still sizeable, at 
the 10% level for large Mi nv . The large effects at high pt and Mi nv are dominated by the 
double logarithmic contributions arising from virtual W and Z bosons in loop diagrams. 

The small peaks visible in the gg invariant mass distribution correspond to two-particle 
thresholds related to b\b\, b^b2, and pairs in gg vertex and box diagrams, illustrated 
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channel 


full result 


PT < lOU KjQ V 


Pt < ^ou Lre v 


(SPS la') 


at NLO [fb] 


& M < 2.5 [fb] 


& M < 2.5 [fb] 


99 


3280 


1643 (-50%) 


778 (-76%) 


qq 


427 


373 (-13%) 


280 (-34%) 


91 


58.5 


30.6 (-48%) 


16.2 (-72%) 



Table 2: Integrated hadronic cross section at NLO within the SPS la' scenario for the different 
production channels. Comparison of the full (unconstrained) results and cross sections where cuts 
on the pseudo-rapidities r\ and on the transverse momenta pt of the outgoing top-squarks are 
applied. The relative changes compared to the full results are given in brackets. 



in Fig. A. 4 [in the SPS la' scenario, the masses of the involved squarks are = 460.7 GeV, 
= 514.8 GeV, m ( - 2 = 569.4 GeV]. Thresholds from the squarks of the first two genera- 
tions are CKM suppressed. The threshold effects appear also in the p^-distribution, around 
300 GeV, but they are smeared out and much less pronounced. 

Fig. | shows total K factors, defined as K = {a^ g LO + a^ LO + a^)/{a^ + <jg>). It is 
obvious that, although small for the total cross section, the EW higher order contributions 
cannot be neglected for differential distributions where, in the high-p^ and high-Mj nt , range, 
they are of the same order of magnitude as the SUSY-QCD corrections (H 




Figure 4: Same as Fig. |^, but shown are the K factors, K = a NLO /a LO , for gg fusion (upper 
plots) and qq channels (lower plots). 




Figure 5: Same as Fig. §, but shown is the total K factor, K = (cr^ LO +cr^ LO +<?g°)/{Vg° +^°)- 



4.3 SUSY parameter dependence 

In order to study the dependence of the EW contributions on the various SUSY parameters 
in more detail, we consider the ratio of the NLO contribution in each channel to the com- 
bined gg + qq Born cross section, 5tot = ^ a {gg qq gi )l a tot ■ We focus on those parameters 
that determine the top-squark mass, cf. Eq. ([O]), and vary each quantity out of the set 
m <5 3 ' m t7 3 ' tan/3, A t , or /u around its SPS la' value while keeping all other parameters 
fixed to those of the default SPS la' scenario. The results are displayed in the left panels 
of Figs.^ - |l^. Simultaneously, we show the mass of the light top-squark i\ as a function 
of the varied parameter in the respective right panels (black solid lines). The parameter 
configuration of the SPS la' scenario is marked by a vertical gray dotted line in all the 
figures. 

We find the following general behaviors. The 57 contributions are from tree level dia- 
grams and the only relevant parameter is thus the top-squark mass m^. In all scenarios, 
the 57 fusion channel is as important as the EW corrections to the qq and gg processes. The 
qq corrections, being practically always negative, involve many different SUSY particles in 
the loops, although the relative corrections show only small variations. The gg contribu- 
tions are more sensitive to the considered SUSY parameters. The plots show striking peaks 
(some of them are also visible in qq annihilation), which correspond to threshold effects and 
can be explained by the SUSY particle masses displayed at the right panels of Figs. [6| - 10 



They occur in the Higgs-exchange diagrams when = m H o/2 (red long-dashed lines in 
the figures), and in the top-squark wave function renormalization when equals the sum 
of masses of a neutralino and the top-quark (green dash-dotted lines) or of a chargino and 
the bottom-quark (blue dashed lines). The chargino- induced peaks are less pronounced 
than those from neutralinos and not visible in Fig. |?j and Fig. 1C . 

Outside of such singular parameter configurations, over a wide range of SUSY pa- 
rameters, the combined EW contributions to top-squark pair production are only weakly 
parameter dependent. 
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Figure 6: Left: Relative EW corrections as a function of the soft-breaking parameter TOq 3 for 
each of the indicated channels compared to the combined (gg + qq) LO cross section in the SPS la' 
scenario where is varied around the SPSla' value (gray dotted line). Right: Mass of t\, half 
of the mass of H°, sums of the masses of the top-quark and Xi> X°j an d X§ , respectively, and sum 
of the masses of the bottom-quark and as a function of ?71q 3 ■ All other parameters are chosen 
according to the SPSla' scenario. 
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Figure 7: Same as Fig. ^, but for variation of the soft-breaking parameter m 
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Figure 8: Same as Fig. [| but for variation of tan/3. 
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Figure 9: Same as Fig. [?], but for variation of trilinear coupling parameter A t 
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Figure 10: Same as Fig. |6[ but for variation of the Higgs parameter 
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5. Conclusions 

We have completed the NLO calculation for the tt* production at hadron colliders by 
providing the complete EW corrections at the one-loop level. 

To obtain a consistent and IR- finite result, we have considered the interference terms 
between QCD and EW NLO terms for both virtual and real contributions. Also, a new 
class of photon-induced partonic processes of tt* production occurs, which was found to 
yield considerable contributions, comparable in size to the corrections to qq annihilation 
and gg fusion or even larger. 

In total, the NLO EW contributions reach in size the 10-20% level in the pr and 
invariant-mass distributions and are thus significant. Outside singular parameter config- 
urations associated with thresholds, the dependence on the MSSM parameters is rather 
smooth. 



Recently, a preprint appeared on the same topic [32], where the authors consider virtual 
corrections and the soft part of the real corrections, both for the gg fusion channel; the 
hard part of the real corrections, as well as the contributions from the other channels are 
missing. The numerical results can therefore not directly be compared with ours at this 
stage. 
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Appendix 

A. Feynman diagrams 

We show here generic Feynman Diagrams for the pair production of lighter top-squark at 
0{aa 2 s ). Diagrams for t^ti production can be constructed in complete analogy. The qq 
annihilation channels are exemplified by uu annihilation. Furthermore, the label S refers 
to all neutral Higgs (and Goldstone) bosons hP, H°, A , G°, and the label S to all charged 
Higgs (and Goldstone) bosons H^, G ± . 
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Figure A.l: Feynman diagrams for top-squark pair production at the Born level via gg fusion 
(left) and qq annihilation (right), here shown for u-quarks. As in the following figures, diagrams 
with crossed final states are not shown explicitely. 




Figure A. 2: Feynman diagrams for gluon-photon fusion. 




Figure A. 3: Feynman diagrams for virtual corrections to top-squark pair production via qq an- 
nihilation (here for u-quarks). The label S° refers to all neutral Higgs bosons h°, H°, A , G°, the 
label S to all charged Higgs bosons H ± , G ± . (a) counter-term diagrams, (b) vertex corrections, 
(c) IR singular box diagrams, (d) IR finite box diagram. 
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Figure A. 4: Fcynman diagrams for virtual corrections to top-squark pair production via gg fusion, 
diagrams with crossed final states are not explicitely shown. The label S° refers to all neutral Higgs 
bosons h°, H°, A , G°, the label 5 to all charged Higgs bosons H ± , GA (a) counter-term diagrams, 
(b) vertex corrections, (c) self-energy corrections, (d) box diagrams. 
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Figure A. 5: Feynman diagrams for real photon radiation, (a) IR divergent - (b) IR finite contribu- 
tions for the gg channel; (c) IR divergent - (d) IR finite contributions for the qq channels. Feynman 
diagrams with photon radiation off the other quark or squark and with crossed final states arc not 
shown explicitely. 




Figure A. 6: Feynman diagrams for gluon bremsstrahlung from the QCD and EW Born diagrams 
(radiation from upper legs is not explicitly shown). Only interference terms between initial and 
final state gluon radiation are non- vanishing. 




Figure A. 7: Feynman diagrams for box contributions of 0(a 2 s ) (left) interfering with electroweak 
Born graphs (right), here for uu annihilation. 



- 23 - 



References 



J. R. Ellis and S. Rudaz, Search for super symmetry in toponium decays, Phys. Lett. B128 
(1983) 248. 



[1] 

[2] A. Djouadi, J. Kalinowski, P. Ohmann and P. M. Zerwas, Heavy SUSY Higgs bosons at e + e~ 
linear colliders, Z. Phys. C74 (1997) 93-111 fhep-ph/9605339|| . 

[3] ALEPH Collaboration, A. Heister et. al, Search for scalar quarks in e + e~ collisions at y/s 



up to 209 GeV, Phys. Lett. B537 (2002) 5-20 |hep-ex/0204036 ]. 

DELPHI Collaboration, P. Abreu et. al., Search for super symmetric partners of top and 



bottom quarks at *fs = 189 GeV, Phys. Lett. B496 (2000) 59-75 jhep-ex/0103034 |. 

L3 Collaboration, M. Acciarri et. al., Searches for scalar quarks in e + e~ interactions at ^fs 



= 189 GeV, Phys. Lett. B471 (1999) 308-320 piep-ex/9910020|l . 

OPAL Collaboration, G. Abbiendi et. al., Search for scalar top and scalar bottom quarks at 



LEP, Phys. Lett. B545 (2002) 272-284 | hep-ex/0209026 ] 



[4] A. C. Kraan, Stop and sbottom searches at LEP, [hep-ex/0305051] 

[5] 



CDF Collaboration, A. A. Affolder et. al., Search for scalar top quark production in pp 
collisions at = 1.8 TeV, Phys. Rev. Lett. 84 (2000) 5273-5278 lhep-ex/991201S| ]. 

DO Collaboration, V. M. Abazov et. al., Search for pair production of light scalar top quarks 
inpp collisions at s/s = 1.8 TeV, Phys. Rev. Lett. 93 (2004) 011801 |hep-ex/0404028(| . 



[6] CDF Collaboration, T. Aaltonen, Search for Direct Pair Production of Supersymmetric Top 



and Supersymmetric Bottom Quarks inpp Collisions at yj\s) — 1.96 TeV, arXiv : 0707 . 2567 



[hep-ex] 



DO Collaboration, V. M. Abazov et. al., Search for the pair production of scalar top quarks in 
the acoplanar charm jet topology inpp collisions at y/s = 1.96-TeV, Phys. Lett. B645 (2007) 
119-127 ]hep-ex/0611003| |. 

DO Collaboration, V. M. Abazov et. al., Search for squarks and gluinos in events with jets 
and missing transverse energy inpp collisions at y/s — 1.96-TeV, Phys. Lett. B638 (2006) 
119-127 jhep-ex/0604029| . 



[7] 
[8] 



CDF and DO Collaboration, T. Nunncmann, Searches for scalar top and bottom quarks at 
the Tevatron, PoS HEP2005 (2006) 348 |iep-ex/0602038| . 



HI Collaboration, A. Aktas et. al., Search for bosonic stop decays in R-parity violating 
supersymmetry in e + p collisions at HERA, Phys. Lett. B599 (2004) 159-172 



[ hep-ex/0405070 1. 

ZEUS Collaboration, S. Chekanov et. al., Search for stop production in R-parity-violating 



[9] 



supersymmetry at HERA, Eur. Phys. J. C50 (2007) 269-281 ftiep-ex/061101S| ]. 

G. L. Kane and J. P. Lcvcillc, Experimental constraints on gluino masses and supersymmetric 
theories, Phys. Lett. B112 (1982) 227. 

P. R. Harrison and C. H. Llewellyn Smith, Hadroproduction of supersymmetric particles, 
Nucl. Phys. B213 (1983) 223. 

E. Reya and D. P. Roy, Supersymmetric particle production at pp collider energies, Phys. 
Rev. D32 (1985) 645. 

S. Dawson, E. Eichten and C. Quigg, Search for supersymmetric particles in hadron - hadron 
collisions, Phys. Rev. D31 (1985) 1581. 



- 24 - 



H. Baer and X. Tata, Component formulae for hadroproduction of left-handed and 
right-handed squarks, Phys. Lett. B160 (1985) 159. 

[10] W. Beenakker, R. Hopker, M. Spira and P. M. Zerwas, Squark production at the Tevatron, 



Phys. Rev. Lett. 74 (1995) 2905-2908 Hhep-ph/9412272 ]. 

W. Beenakker, R. Hopker, M. Spira and P. M. Zerwas, Squark and gluino production at 



hadron colliders, Nucl. Phys. B492 (1997) 51-103 [tiep-ph/9610490 



[11] 
[12] 
[13] 
[14] 
[15] 



W. Beenakker, M. Kramer, T. Plehn, M. Spira and P. M. Zerwas, Stop production at hadron 
colliders, Nucl. Phys. B515 (1998) 3-14 |hep-ph/971045l|| . 



G. Bozzi, B. Fuks and M. Klasen, Non-diagonal and mixed squark production at hadron 



colliders, Phys. Rev. D72 (2005) 035016 |hep-ph/0507073| |. 

A. Bartl et. ai, Search of stop, sbottom, tau-sneutrino, and stau at an e + e~ linear collider 



with = 0.5-TeV to 2 TeV, Z. Phys. C76 (1997) 549-560 |hep-ph/9701336 



H. E. Haber and G. L. Kane, The Search for Super symmetry: Probing Physics Beyond the 
Standard Model, Phys. Rept. 117 (1985) 75-263. 

J. Kublbeck, M. Bohm and A. Denner, FeynArts: Computer Algebraic Generation of 
Feynman Graphs and Amplitudes, Comput. Phys. Commun. 60 (1990) 165-180. 

T. Hahn, Generating Feynman diagrams and amplitudes with FeynArts 3, Comput. Phys. 
Commun. 140 (2001) 418-431 friep-ph/0012260 . 



T. Hahn and C. Schappacher, The implementation of the minimal supersymmetric standard 
model in FeynArts and FormCalc, Comput. Phys. Commun. 143 (2002) 54-68 



[hep-ph/0105349 



[16] T. Hahn and M. Perez- Victoria, Automatized one-loop calculations in four and D dimensions, 



Comput. Phys. Commun. 118 (1999) 153-165 |)iep-ph/9807565 |. 

T. Hahn and M. Rauch, News from FormCalc and LoopTools, Nucl. Phys. Proc. Suppl. 157 



(2006) 236-240 |hep-pri/0601248[ 



[17] 



[18] 



[19] 



G. 't Hooft and M. J. G. Veltman, Scalar One Loop Integrals, Nucl. Phys. B153 (1979) 
365-401. 

G. Passarino and M. J. G. Veltman, One loop corrections for e+ e- annihilation into mu+ 
mu- in the Weinberg model, Nucl. Phys. B160 (1979) 151. 

W. Beenakker and A. Denner, Infrared Divergent Scalar Box Integrals With Applications in 
the Electroweak Standard Model, Nucl. Phys. B338 (1990) 349-370. 

A. Denner, U. Nierste and R. Scharf, A Compact expression for the scalar one loop four point 
function, Nucl. Phys. B367 (1991) 637-656. 

S. Heinemeyer, W. Hollik and G. Weiglein, FeynHiggs: A program for the calculation of the 
masses of the neutral CP-even Higgs bosons in the MSSM, Comput. Phys. Commun. 124 



(2000) 76-89 jhep-ph/9812320 



T. Hahn et. al., Higher-order corrected Higgs bosons in FeynHiggs 2.5, [hep-ph/0611373] 



F. Bloch and A. Nordsieck, Note on the radiation field of the electron, Phys. Rev. 52 (1937) 
54-59. 



[20] 

[21] W. Beenakker, S. C. van der Marck and W. Hollik, e + e~ annihilation into heavy fermion 



pairs at high-energy colliders, Nucl. Phys. B365 (1991) 24-78. 



- 25 - 



U. Baur, S. Keller and D. Wackeroth, Electroweak radiative corrections to W boson 



production in hadronic collisions, Phys. Rev. D59 (1999) 013002 hep-ph/ 9807417 1 . 

S. Dittmaier and M. Kramer, Electroweak radiative corrections to W-boson production at 
hadron colliders, Phys. Rev. D65 (2002) 073007 friep-ph/0109062| ]. 

G. Altarelli and G. Parisi, Asymptotic Freedom in Parton Language, Nucl. Phys. B126 
(1977) 298. 

D. Wackeroth and W. Hollik, Electroweak radiative corrections to resonant charged gauge 
boson production, Phys. Rev. D55 (1997) 6788-6818 [ frep-ph/9606398| ]. 



K. P. O. Diener, S. Dittmaier and W. Hollik, Electroweak radiative corrections to 
deep-inelastic neutrino scattering: Implications for NuTeV?, Phys. Rev. D69 (2004) 073005 



[ hep-ph/0310364 ]. 

A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Parton distributions 



incorporating QED contributions, Eur. Phys. J. C39 (2005) 155-161 |hep-ph/0411040 |. 
O. Brein and W. Hollik, Distributions for MSSM Higgs boson + jet production at hadron 



colliders, Phys. Rev. D76 (2007) 035002 [arXiv : 0705 . 2744 [hep-ph] 



J. A. Aguilar-Saavedra et. ai, Supersymmetry parameter analysis: SPA convention and 



project, Eur. Phys. J. C46 (2006) 43-60 |hep-ph/0511344 l. 
See: http://www-cdf.fnal.gov/physics/new/top/top.html. 

B. C. Allanach et. at, The Snowmass points and slopes: Benchmarks for SUSY searches, 



[hep-ph/0202233] 



[32] 



http://www.ippp.dur.ac.uk/~georg/sps/sps.html. 

M. Beccaria, G. Macorini, L. Panizzi, F. M. Renard and C. Verzegnassi, A survey of one-loop 



electroweak super symmetric effects in stop-antistop production at LHC, arXiv : 0710 . 5357 
[hep-ph] 



- 26 - 



